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Heterogeneous adoption thresholds exist widely in social contagions, but were always neglected in previous 
studies. We first propose a non-Markovian spreading threshold model with general adoption threshold distri¬ 
bution. In order to understand the effects of heterogeneous adoption thresholds quantitatively, an edge-based 
compartmental theory is developed for the proposed model. We use a binary spreading threshold model as a 
specific example, in which some individuals have a low adoption threshold (i.e., activists) while the remaining 
ones hold a relatively high adoption threshold (i.e., bigots), to demonstrate that heterogeneous adoption thresh¬ 
olds markedly affect the final adoption size and phase transition. Interestingly, the first-order, second-order and 
hybrid phase transitions can be found in the system. More importantly, there are two different kinds of crossover 
phenomena in phase transition for distinct values of bigots’ adoption threshold: a change from first-order or hy¬ 
brid phase transition to the second-order phase transition. The theoretical predictions based on the suggested 
theory agree very well with the results of numerical simulations. 
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I. INTRODUCTION 

Social contagions are ubiquitous in human society, and gen¬ 
erate new scientific challenges for network science iniEi. 
All the researches about sentiments contagion, information 
spreading and behavior spreading fall into the category of 
studying social contagions IS a. In particular, behavior 
spreading, as a representative and essential type of social con¬ 
tagions, has attracted great attention, both theoretically and 
experimentally Ea. Understanding the spreading mecha¬ 
nisms behind behavior has the potential to not only help us 
design better anti-virus strategies ||71, but also shed new in¬ 
sights into the control of social unrest |[T]. Moreover, different 
from biological contagions (such as epidemic spreading) Ha- 
Ha, social contagions display one inherent characteristic: so¬ 
cial reinforcement effect E fT4ll . which usually plays a vital 
role in the final status of contagions. 

To examine social reinforcement effect, some successful 
models have been proposed iiiaiia, most of which assume 
that all individuals have the same adoption threshold to incor¬ 
porate such an effect (so-called threshold model). That is to 
say, each individual adopts the behavior only when the frac¬ 
tion oa or number UB of neighbors with the adopted state 
exceeds his adoption threshold. In this case, the social conta¬ 
gions is a trivial case of Markovian process. By means of nu¬ 
merical simulations and theoretical analysis, it was found that 
the social reinforcement effect can evidently change the phase 
transitions of contagion dynamics. More specifically, the fi¬ 
nal adoption size first grows continuously and then decreases 
discontinuously with the increase of mean degree when the 
adoption threshold of all individuals is identical Ca. Af¬ 
ter this seminal discovery, the role of threshold model and 
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its various underlying mechanisms, in social contagions have 
been intensively explored, including the influence of dynami¬ 
cal parameters (e.g., initial seeds and threshold sizes iniiia), 
and topology characteristics (degree-degree correlations QSl, 
clustering MM, community structure m as well as mul- 
tiplexity framework ll23l l24l ). which are the primary factors 
in determining the final adoption size and phase transition. 
Some non-Markovian social contagion models were also pro¬ 
posed to describe the social reinforcement effect Il25lj27ll . For 
example, in a recent research paper ll26l . where the social re¬ 
inforcement was derived from the memory of non-redundant 
information transmission, it was found that the perturbation 
of dynamical or structural parameters makes the dependence 
of final behavior adoption size on information transmission 
probability change from being discontinuous to being contin¬ 
uous. In Ref. EH, Zheng et al. further verified that the role of 
social reinforcement in behavior diffusion on regular graphs 
and online social networks, which is consistent with early ex¬ 
perimental anticipation ll2^ 

Recently, researchers found that the widely existed indi¬ 
vidual heterogeneity dramatically alters the spreading dynam¬ 
ics. In epidemic spreading, the heterogeneous infectivity 
and susceptibility can change the outbreak threshold ll29l[MI . 
For information spreading, the heterogeneous waiting and re¬ 
sponse time may speed up or slow down the velocity of in¬ 
formation diffusion ll3n[32]| . Statistical physicists found that 
individual heterogeneity can induce the hybrid phase tran¬ 
sition lf3^ lT4l . which mixes the traditional first-order and 
second-order transitions, in k-core percolation 051 and boot¬ 
strap percolation IITtIO^OtII . In practical behavior spreading, 
individuals usually show different wills to mimic the behavior, 
which means that each agent owns his own adoption thresh¬ 
old OSl . Some individuals with low adoption threshold show 
strong wills to adopt the behavior and act as activists. Nev¬ 
ertheless, others with high adoption threshold need to capture 
more behavioral information before imitation, and they often 
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act as bigots. With regard to the difference of adoption thresh¬ 
old, it may be closely related with personal interests, educa¬ 
tion background, or other personality and social factors HI. 
For example, the well-educated population is more likely to 
adopt the high-tech products than that among the populations 
who lacks the basic education. Similarly, students are more 
likely to adopt an interesting computer game than housewives. 

Unfortunately, there is still absence of systematical under¬ 
standing about the role of heterogeneous adoption thresholds 
in social contagions. Aiming to resolve this issue, we will ex¬ 
plore how heterogeneous adoption thresholds affect the final 
adoption size and phase transition of social contagions based 
on a so-called binary spreading threshold model, which is a 
non-Markovian process. Meanwhile, an edge-based compart- 
mental theory is developed for quantitative validation. Inter¬ 
estingly, it is found that heterogeneous adoption thresholds 
significantly affect the final adoption size, and generate the 
hierarchical characteristic of adopting behavior; activists first 
adopt the given behavior themselves, and then stimulate bigots 
to follow this behavior. Moreover, it is worth noting that such 
a heterogeneous threshold model results in the existence of 
first-order, seconde-order and hybrid phase transitions. More 
importantly, there are two different kinds of crossover phe¬ 
nomena in phase transition for distinct values of the bigots’ 
adoption threshold; a change from first-order or hybrid phase 
transitions to the second-order phase transition. In what fol¬ 
lows, we will first describe the heterogeneous social conta¬ 
gion model in complex networks, followed by the description 
of edge-based compartmental theory, and then represent the 
simulation and analysis results. Finally, we will draw our con¬ 
clusions. 


II. SOCIAL CONTAGION MODEL WITH 
HETEROGENEOUS ADOPTION THRESHOLDS 

Behavior spreading on complex networks is considered 
with N nodes and the degree distribution P{k). As the in¬ 
teraction networks, we use the configuration model ll39l to 
avoid the additional influence of degree-degree correlations. 
Nodes in the network represent individuals and edges between 
nodes stand for the contacts with which behavioral informa¬ 
tion transmission may occur. For each individual, a static be¬ 
havioral adoption threshold is assigned according to a spe¬ 
cific distribution function F{T), which is independent of net¬ 
work topology. The larger value of T means that an individual 
needs to capture more behavioral information from his neigh¬ 
bors before adopting the behavior. 

With regard to behavior spreading dynamics, we generalize 
the spreading threshold model with social reinforcement de¬ 
rived from memory of non-redundant information transmis¬ 
sion characters ll26ll40l . In this model, each individual falls 
into one of the three states; susceptible (S), adopted (A) and 
recovered (R) (namely, susceptible-adopted-recovered, SAR 
model). In the susceptible state, an individual does not adopt 
the behavior. In the adopted state, an individual adopts the be¬ 
havior and tries to transmit the behavioral information to his 
neighbors. In the recovered state, an individual loses interest 


in the behavior and will not transmit the behavioral informa¬ 
tion further. Initially, a vanishingly small fraction of individ¬ 
uals pq are chosen as seeds (adopters) at random, while the 
others are fixed in the susceptible state. 

At each time step, each adopted individual v tries to dif¬ 
fuse the behavioral information to every susceptible neigh¬ 
bor with probability A. In particular, once the information 
is transmitted through an edge successfully, it will never be 
transmitted again, i.e., only non-redundant information trans¬ 
mission is allowed. If the susceptible neighbor u of u is suc¬ 
cessfully informed, his cumulative pieces of information m 
add 1 (i.e., m — m + 1). Subsequently, individual u com¬ 
pares the new value of m with his adoption threshold T^, and 
becomes an adopter once m>Tu. Obviously, whether an in¬ 
dividual adopts the behavior is determined by the cumulative 
pieces of behavioral information he ever received from dis¬ 
tinct neighbors. Thus, the non-Markovian effect is induced in 
the behavior spreading dynamics. After information transmis¬ 
sion, individual v may lose interest in the behavior with prob¬ 
ability 7 and then moves into the recovered state. Individuals 
falling into the recovered state will stop from participating in 
the further behavioral information spreading, and the spread¬ 
ing dynamics terminate when all adopted individuals become 
recovered. 


III. EDGE-BASED COMPARTMENTAL THEORY 

In order to describe the strong dynamical correlations 
among the states of neighbors in heterogeneous social con¬ 
tagion model, an edge-based compartmental approach is es¬ 
tablished herein, which is inspired by Refs. EJil- Cor¬ 
respondingly, the notations S{t), A(t) and R(t) respec¬ 
tively represent the fraction of individuals in the susceptible, 
adopted, and recovered states at time step t. 


A. General adoption threshold distrihution 

Individual u is set to be in the cavity state, which means that 
he can receive behavioral information from neighbors but not 
transmit behavioral information to his neighbors Il44l . Define 
6{t) as the probability that individual v has not transmitted the 
behavioral information to individual u along a randomly cho¬ 
sen edge by time t. Thus, the probability that an individual u 
with degree k has received m pieces of behavioral information 
from distinct neighbors by time t can be expressed as 

= ( 1 ) 

Individual u in the susceptible state implies that the cumula¬ 
tive pieces of behavioral information m he received are still 
less than his adoption threshold T„. According to the social 
contagion model in Sec.[^ the adoption threshold and degree 
are independent. Considering all possible values of m and 
T„, it can be obtained that the probability of individual u with 
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degree k being susceptible is 


Combining Eqs. (j^ and (j^, one has that 


T„-l 

s{k,t) = '^F{Tu) (l)m.{k,t). ( 2 ) 

Tu m—0 

Combing the degree distribution of a network, the fraction of 
susceptible individuals at time step t is 

S{t)=J2Pik)s{k,t). (3) 

k 

It is obvious that S{t) can be figured out after 0{t) is known. 

As a neighbor of individual u may be in one of susceptible, 
adopted or recovered state, 9{t) can be divided into three cases 
as 


(^{t) = is{t) + ^Ait) + ^R{t). (4) 

And ^s(f) [CA(i) or denotes that the susceptible 

(adopted or recovered) neighbor u of m has not transmitted 
the behavioral information to individual u up to time t. 

Then, let’s explore the above three terms. If individual v 
with degree k' is susceptible initially, he can not transmit the 
behavioral information to u, but only receive it from other k' — 
1 neighbors since u in a cavity state. Thus, the probability that 
individual v has received m pieces of behavioral information 
by time t is 

T^{k',t) = (^' (5) 

Taking all possible values of m and T„ into consideration, the 
probability of individual v remaining in susceptible state be¬ 
comes [similar to Eq. 0 ] 

T„-l 

e{k',t) = ^ T^{k',t). ( 6 ) 

T„ m=0 


= 


7 [l- 6 i(f)](l-A) 
A 


( 10 ) 


Inserting Eqs. (j7]l and (10 1 into Eq. Q, the following expres¬ 
sion can be obtained that 


Uit) = 0{t)- 


Ek'k'Pik')Q(k',t) 7[l-0(t)](l-A) 


ik) 


A 


( 11 ) 


Substituting Eq. (Ill into Eq. ([^, we get the time evolution of 
6(t) in details 




dt 


(fc) 


( 12 ) 

Susceptible individuals move into the adopted states once 
they adopt the behavior, meanwhile the adopted individuals 
may lose interest in the behavior and become recovered. Thus, 
we can easily get the evolution of adopted and recovered indi¬ 
viduals as 


dA{t) dS{t) 


dt 


dt 


- iMt) 


(13) 


and 


dR{t) 

dt 


= iMt), 


(14) 


respectively. Eqs. ([2l-([^ and ([T^-([l4)l give us a complete and 
general description of heterogeneous social contagions, from 
which the fraction in each state at arbitrary time step can be 
calculated. When f —c», we can get the final adoption size 

R{oo). 


B. Binominal adoption threshold distribution 


The probability of an edge connecting an individual with de¬ 
gree k' is k'P{k')/(k) for uncorrelated networks, where (k) 
is the mean degree. Thus, it can be obtained that 


^ k'P{k')Q{k'A) 


(7) 


Subsequently, we turn to the expressions of and 

^R(t). Once the behavioral information is transmitted through 
an edge with probability A, the edge will no longer satisfy the 
definition of 0{t). Thus, the rate of flow from 6{t) to 1 — 9(t) 
is which can be expressed as 


d9{t) 

dt 


-X^Ai.t). 


( 8 ) 


Eor the growing of ^R{t), two events conditions must be met 
simultaneously. At time t, the behavioral information does 
not transmit through an edge with probability 1 — A, and the 
adopted individual enters recovered state with probability 7 . 
Then, 


d^R{t) 


In this subsection, we pay attention to the behavior adoption 
threshold with a binominal distribution F{T). More specifi¬ 
cally, a fraction of individuals p have a relatively low adoption 
threshold Ta, whereas the remaining individuals have a high 
adoption threshold Tb. F{T) can be expressed as 


F{T) = 


Ta, with probability p, 
Tb, with probability 1 — p. 


(15) 


Eor simplicity, the values of adoption thresholds are defined as 
Ta = 1 and Tb > 1. Individuals with low adoption threshold 
Ta are considered as activists, while those with high adop¬ 
tion threshold Tb are regarded as bigots. We herein name this 
kind of social contagion model as binary spreading threshold 
model. And the edge-based compartmental theory is utilized 
to analyze the binary spreading threshold model by substitut¬ 
ing Eq. into various equations that give the solutions of 
S{t), A{t) and R{t). Particularly, we rewrite Eqs. (j^ and (j^ 
as 


n-i 

s{k,t) =p9{tf -f (1 -p) ^ (j)m{k,t), (16) 

m=0 


dt 


7(1 - 


(9) 
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Q{k\t)=pe{tf ^ + (1-p) ^ (17) 

m—0 

respectively. At time t, the fractions of susceptible individuals 
in the activist and bigot populations are given by 

k 

and 



6(oo) 0(°°) 


n-i 

Sh{t)=Y,Pik)Y,^ra{k,t), (19) 

k m—0 

respectively. Considering the fractions of the activist and 
bigot populations, the density of susceptible individuals at 
time step t can also be written as 

S{t) = pSi{t) + {I - p)Sh{t). (20) 


The effects of heterogeneous adoption threshold on phase 
transition is another issue concerned. To analyze the phase 
transition, we can address the fixed point (root) of Eq. © at 
the steady state (i.e., t oo) with Eq. (17 1 . That is the fixed 
point of 


0(oo) = J/[0((X))], (21) 

where 

7[1-6»(c3o)](1-A) 

( 22 ) 

Erom Eig. [Ja), it can be seen that the number of nontrivial 
roots is either 0, 1 or 3 when > 3. With 7), = 2, the 
number of nontrivial roots is either 0 , 1 or 2 [see Eig.[T][b)]. 


1. Case ofTb > 3 


In this subsection, we discuss the case of > 3. Eor the 


given P{k), p and 7 , Eq. (21 1 has only one trivial solution 


0{oo) = 1 when A is small. With the increase of A, 9{oo) 
decreases continuously to a nontrivial solution 9{oo) < 1 first 
[see example in Eig. [^a)], which means that R{oo) grows 
continuously first. That is to say, there is a second-order (con¬ 
tinuous) phase transition. By setting 9{oo) and y[0(c»)] tan¬ 
gent at 0 ( 00 ) = 1 lIMl l45l . we get the continuous critical 
information transmission probability as 


xi' = 


l{k) 


p((fc2)-(fc))-(l-7)(fc)^ 


(23) 


where (fc) and {k"^) are the first and second moments of degree 
distribution, respectively. The critical value separates the 
local behavior adoption (i.e., the behavior can be adopted by 
a vanishingly small fraction of individuals) from the global 


FIG. 1. (Color online) Illustration of graphical solutions of Eq. ID 
when Tj, = 4 (a) and Tj, = 2 (b) on random regular networks. In 
(a) and (b), the black solid lines represent the horizontal axis and the 
red and green dots denote the tangent points. The inset of (a) and (b) 
shows the physically meaningful solutions of Eq. 0 when Ti, = 4 
and Tb = 2, respectively. And the arrows indicate critical transmis¬ 
sion probabilities. The fixed points of Eq. l |21^ are the intersections 
between the curve and horizontal axis. Other parameters are defined 
p — 0.3 and (fc) = 10. 


behavior adoption (i.e., the behavior can be adopted by a fi¬ 
nite fraction of individuals). Erom Eq. (j2^, it is discovered 
that the occurrence of global behavior adoption is determined 
by the network topology (i.e., degree distribution), the frac¬ 
tion of activists p and the recovery probability 7 . The global 
behavior adoption is more likely outbreak (i.e., a lower X^J) 
for scale-free networks with divergent second moment degree 
distribution [i.e., (fc^) —1 00 ]. Increasing the value of p can fa¬ 
cilitate the global behavior adoption (i.e., a lower X^J). When 
p = 1, Eq. ( [23| returns to the case of epidemic outbreak 
threshold B^ . 


Eixing all the parameters except p, similar to Eq. (23 1 , we 
get the continuous critical fraction of activists as 


Pi^ = 


(1 - 7 )( fc ) 

X{{B)-{k)) ' {k^)-{k)- 


l{k) 


+ 


(24) 


Erom Eq. it can be known that an enough fraction of ac¬ 
tivists are necessary for triggering the global behavior adop¬ 
tion, and p^J decreases with the increase of network hetero¬ 


geneity, p and A. By setting 7 = 0 in Eq. (24 1 , another crit¬ 
ical proportion of activists p* can be obtained, below which 
any values of A can not trigger the global behavior adoption, 
which is given by 


P 


★ 

C 


(fc) 

(P)-(fc)- 


(25) 


It is worth noting that Eq. ( [25] ) is the same with network’s 
percolation condition na, which means that the global be¬ 
havior adoption is possible only when the activists percolate 
the entire network (i.e., activists can form a finite connected 
cluster). 

As shown in Eig. [TJa), three nontrivial roots of Eq. (|2T| oc¬ 
cur when A is large enough [see Eig. [^a) for A = 0.615]. 
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This phenomenon is caused by the bigots, since Eq. ( [2T] i has 
at most one nontrivial root when the bigots is absent |26|. In 
this case, the meaningful solution will be given by the largest 
stable root (since only this value can be achieved physically). 
For A = Ag = 0.635, the tangent point is the solution. For 
A > Ag, the meaningful solution is the only stable fixed 
point. The meaning solution of Eq. (j2T| changes abruptly to a 
small value from a relatively large value at A^ [see the insert 
Fig. [^a)], and resulting in a discontinuous growth of R{oo). 
Based on bifurcation theory ll46l . the discontinuous critical 
information transmission probability is gained as follows 


I ^ 7 

A + 7 - 1 ’ 


(26) 


where 


A = 


{k) 


and 9s{oo) is the fixed point of Eq. (21 1 . Combining Eqs. (|^ 
and ( [T7 ] i, we get 

dQ{k', oo) 


= p(fc'-l)0(oo)'='-2 + (l-p)A’, (27) 


d0{t) 


where 


1 /If _ 1 \ 

T' = ^ ( 1 {{k' -m- 1 ) 0 ( 00 )'= - 0 ( 00 )]"* 

m=0 A / 

- me{oof-^-^[l - 0(oo)]'"-i}. 

(28) 


Using the analytical method similar to Eq. ( [2^ , the discontin¬ 
uous critical fraction of activists can be expressed as 


7, A(l-y)+7(l-A) 
A(Z - j;) 


(29) 


2. Case ofTb = 2 


We study the special case of T;, = 2 in this subsection. 
As shown in Fig. for any A, y[0(oo)] can only tan¬ 
gent to 0(oo) when 0(oo) = 1, and can not tangent to other 
0(oo) < 1. In addition, Eq. (21 1 has 1 or 2 nontrivial fixed 
points. These phenomena means that the meaningful solution 
of Eq. ( [2 T] i jumps to another one at the critical information 
transmission probability [see the inset of Fig. [TJb)]. As a re¬ 
sult, R{oo) increases discontinuously versus A. The critical 
information transmission probability can be acquired in the 
similar way as when = 3. The values of A^^ and A^ can 
be obtained from Eq. (231, since j/[0(oo)] can only tangent to 
0(oo) when 0(oo) = 1. Numerical solving Eqs. (j2T|l, ( p^ 
and ( [30| , we can get the condition under which the first-order 
phase transition changes to a second-order phase transition. 


IV. NUMERICAL VERIFICATION 


In the study, extensive simulations are conducted for on 
binary spreading threshold model on uncorrelated networks. 
Unless otherwise specified, the network size, mean degree 
and recovery probability are of A = 10,000, {k) = 10 and 
7 = 1.0, respectively. At least 2 x 10^ independent dynam¬ 
ical realizations on a fixed network are used to calculate the 
pertinent average values, which are further averaged over 100 
network realizations. 

The relative variance vr is applied to numerically deter¬ 
mine the size-dependent critical values, such as, A^, A^^, 
and . The relative variance of R{oo) is defined as 


_ {Rjoo) - {R{oo))y 
^ {Rioo)r 


(31) 


where 

Y.u'k'P{k')x 

{k) 

and 


k'{k' - l)P(fc')0'='-2 

¥) 


From the above analysis, we find that R{oo) versus A or 
p first grows continuously and then follows a discontinuous 
fashion. And the continuous and discontinuous growthes of 
R{oo) are caused by the activists and bigots, respectively, 
which can be regarded as hybrid phase transition, mixing the 
traditional first-order and second-order transitions, from the 
perspective of statistical physics nil. Note that the hybrid 
phase transition can change to a second-order phase transi¬ 
tion, and does not exist under any conditions. Numerically 
solving Eqs. (0, (@, and the following equation 


d^g[0(oo)] 

d02(oo) 


(30) 


we can learn the condition under which the hybrid phase tran¬ 
sition disappears. 


where (• • •) denotes ensemble averaging. The value of vr ex¬ 
hibits peaks at the phase transition, which announce the phase 
transition 1481 . We determine the critical value Wc, which rep¬ 
resents Ag, , pI or p^J, as the value of A or p when the 
relative variance reaches its maximum, and 

ujc = argjmax u^}. (32) 

Note that Eq. ( (32) l is not the only possible way to compute the 
critical values, other methods can be used to determine uJc, 
such as, susceptibility Ii49ll50l and variability ED. 


A. Random regular networks 

To be illustrative, we first focus on random regular net¬ 
works (RRNs). Fig.|^a) shows the time evolution of the frac¬ 
tion of adopted individuals Aa{t) and A^jf) in the activist and 
bigot populations with different behavioral information trans¬ 
mission probabilities A. It is found that a hierarchical char¬ 
acter of behavior adoption is caused by heterogeneous adop¬ 
tion thresholds. That is to say, activists with low first 
adopt the behavior and then stimulate the bigots with Ti, to 
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FIG. 2. (Color online) Binary spreading threshold model on random 
regular networks, (a) Average density of adopted individuals Aa{t) 
and Ai,{t) versus time t in the activist and bigot populations with dif¬ 
ferent information transmission probability A. The error bars indicate 
the standard deviations, (b) Final adoption size R{oo) and (c) rela¬ 
tive variance vr versus A for full, activist and bigot populations. In 
figures (a) and (b), symbols represent the simulated results, and the 
lines are the corresponding theoretical predictions from Eqs. un¬ 
til and t^-@. Lines in (c) are the simulation results of vr. In 
(b) and (c), the vertical lines are the critical continuous and discon¬ 
tinuous information transmission probabilities, which are denoted as 
Ac^ and A^, respectively. The theoretical values of A^^ and A^ can be 
gotten from Eqs. and \26\ , respectively. The theoretical (numeri¬ 
cal) values of and A^ separate (b) [(c)] into three regions. Region 
I, only a vanishingly small fraction of individuals can be exposed to 
adopt the behavior; region II, only a finite fraction of activists adopt 
the behavior; region III, a finite fraction of activists and bigots adopt 
the behavior. Other parameters are set to be = 10,000, p = 0.3 
and Tj, = 4. 


adopt the behavior. With a relatively small A = 0.5, Aa{t) 
shows a small peak, and can not stimulate a finite fraction 
of bigots to adopt the behavior [Ai,{t) does not see an ob¬ 
vious peak]. With a relatively large A = 0.7, Aa{t) shows 
a large peak, and further causing the emergence of a large 
peak for Aa{t). The heterogeneous adoption threshold distri¬ 
bution may be used to explain the existence of multimodal in 
the adoption of serves Il38l . The time evolution can be well 
predicted by our edge-based compartmental theory. 

From Figs.|^b) and (c), it can be seen that R{oo) versus 
A shows a hybrid phase transition, which means that R{oo) 
first grows continuously and then follows a discontinuous 
fashion. The continuous and discontinuous phase transitions 
are caused by activists and bigots, respectively. Similar to 


Ref. Il26ll . the discontinuous growth of R{oo) is caused by 
those bigots in the subcritical state who adopt the behavior 
simultaneously. An individual in such a state has received the 
behavioral information but has not yet adopted the behavior, 
and the number of information pieces from distinct neighbors 
is precisely one less than his adoption threshold. The theoret¬ 
ical (numerical) values of A^^ and A^ separate Figs.Sb) [(c)] 
into three regions. The theoretical values of A^^ and can be 
gotten from Eqs. ( [23] l and p6| ), respectively. In region I, with 
A < , both activist and bigot populations adopt the behav¬ 

ior locally (i.e., only a vanishingly small fraction of individu¬ 
als adopted the behavior). In region II, with A^^ < A < A^, 
activists adopt the behavior globally (i.e., a finite fraction of 
activists adopted the behavior) and bigots adopt the behav¬ 
ior locally. In region III with A > A^, both the activists and 
bigots adopt the behavior globally. The numerical values of 
Ag^ and A^ can be obtained by observing vp in Fig.|^c). For 
instance, vp has two peaks, which means that two phase tran¬ 
sitions occur ll48l . And the first peak appears at A^^, while the 
second peak locates at A^. Note that the first (second) peak 
of vp shares the same location with the maximal peak of ac¬ 
tivist population (bigots population <)■ In all, X[, A" 
and R{oo) can be well predicted by our edge-based compart¬ 
mental theory. 


As shown in Fig.[^ R{oo) and the phase transition are sig¬ 
nificantly influenced by the adoption threshold of bigots Tf,. 
And R{oo) decreases with Tb, since a larger value of Tf, re¬ 
quires more information to be exposed for bigots. The phase 
transition is continuous when = 1. In this case, the conta¬ 
gion dynamics is the same with epidemic spreading B5l . The 
phase transition is also continuous when Tt > 6, since there 
are not enough activists to persuade bigots to adopt the be¬ 
havior simultaneously. For the case of Tb = 2, R{oo) shows a 
first-order phase transition, because the bigots are likely to en¬ 
ter subcritical states and adopt the behavior simultaneously. A 
hybrid phase transition emerges with other values of Tb (i.e., 
2 < Tb < 6). As discussed in Sec. IIIB| the type of phase 
transition is verified by bifurcation analysis of Eq. pT| ). In 
Eig. I^b), we verify the phase transition by studying vp in 
simulations. And the simulated values of A^^ and A^ are lo¬ 
cated by studying vp versus A. Eor the second-order phase 
transition, vp has only one peak [see Tb — 1 and Tb — 6 
in Eig. I^c)]. Similarly, vp also has only one peak for the 
first-order phase transition [see = 2 in Eig.|^c)]. Eor the 
hybrid phase transition, vp has two peaks [see 3 < Tf, < 5 in 
Eig.j^c)]. Our theoretical predictions of i?(oo) are well agree 
with the simulation results, except for the cases near the crit¬ 
ical information transmission probability. The deviations be¬ 
tween our predictions and simulation are mainly derived from 
the finite-size effects of the networks, as shown in Eig.|^ The 
deviations of R{oo), A^ and X^J between the simulated and 
theoretical results decrease with network size N. 


Then, we observe R{oo) versus p for different A in Eig. 
We And that R{oo) increases with p. By bifurcation analysis 
of Eq. (21 1 and studying vp, a continuous growth of R{oo) 
is observed with a relatively small A (e.g., A = 0.5), and the 
hybrid phase transition occurs with a relatively large A (e.g., 
A = 0.7 and 0.8). Again, the theoretical and numerical results 
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FIG. 3. (Color online) Effects of bigots’ adoption threshold on bi¬ 
nary spreading threshold model, (a) R(oo) and (b) vr versus A with 
different Tf,. In figure (a), symbols represent the simulated results, 
and the lines are the theoretical predictions from Eqs. and 

-(|201. In figure (b), lines are the simulation results of vr, and we 
plot the value of vr/2Q for Tb = 2. Other parameters are defined as 
N = 10,000 andp = 0.3. 


EIG. 5. (Color online) Effects of the fraction of activists on binary 
spreading threshold model, (a) R{oo) and (b) vr as a function of p 
for different A. In figure (a), symbols represent the simulated results, 
and the lines are the theoretical predictions from Eqs. l|12^-l|14|l and 
cn -(|20^. Lines in (b) are the simulation results of vr. We set other 
parameters to be A'" = 10, 000 and Tb = 4. 
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FIG. 4. (Color online) Finite-size effects on the binary spreading 
threshold model for p = 0.3. (a) R{oo) and (c) vr versus A for 
To = 2 on different network size N. (b) R{oo) and (d) vr versus 
A for Tj, = 4 under different N. The lines in (a) and (b) are the 
theoretical predictions from Eqs. (|12^-l|14[( and (|16)-(|20[(. The lines 
in (c) and (d) are the simulation results of vr. 


agree well. 

From the above analysis, it can be obtained that both A and 
p markedly affect R{oo) and phase transition. Thus, we fur¬ 
ther investigate R{oo) and phase transition on parameter plane 
(A,p) when Tf, = 4 in Fig. Obviously, R{oo) increases 
with A and p. According to the type of phase transition, the 
parameter plane (A,p) is divided into four different regions 
by three vertical lines. The first vertical line can be gotten 
from Eq. ( |25ll, and the other two can be predicted by solving 
Eqs. (21 1 , \26) and (30i. In region I (p < pj = 1/9), there 
are a few activists, who can not percolate the entire popula¬ 
tion. Thus, no matter what the value of A, activists can not be 
made to adopt the behavior globally. When p > p/, the global 
behavior adoption becomes possible, and a crossover phe¬ 
nomenon, which means that the phase transition changes from 
being hybrid to being second-order, occurs in the phase tran¬ 
sition. Meanwhile, the local and global behavior adoptions 
are separated by the red solid curve (i.e., X^J). In region II 
(1/9 < p < 0.15), the relatively few activists lead to the con¬ 
tinuous phase transition. In this region, R{oo) grows continu¬ 
ously versus A for a given p, and a finite fraction of individuals 
adopt the behavior above A/^. With the increase of p, in region 
III (0.15 < p < 0.5), the hybrid phase transition occurs, i.e., 
R{oo) first grows continuously with A and then follows by a 
discontinuous pattern. A finite fraction of activists adopt the 
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FIG. 6. (Color online) Dependence of R{oo) on p and A on ran¬ 
dom regular networks. Color-coded values R{oo) are obtained from 
numerical simulations (a) and theoretical solutions (b) in the param¬ 
eter plane (p, A), and the theoretical predictions are the solutions of 
Eqs. l|12|l-l|14|l and The plane are divided into four re¬ 

gions by three vertical lines, of which the first line can be gotten from 
Eq. 1 25 I and the other two are predicted by solving Eqs. 
and 1 30 '. In region I, only a vanishingly small fraction of individuals 
can be exposed to adopt the behavior (i.e., local behavior adoption). 
Both the regions II and IV show a continuous phase transition, while 
region III exhibits a hybrid phase transition. The red circles (red solid 
curve) and white squares (black dished curve) are the continuous and 
discontinuous simulated (theoretical) critical information transmis¬ 
sion probability, respectively. Moreover, other parameters defined as 
N = 10, 000, p = 0.3 and n = 4. 


FIG. 7. (Color online) Dependence of R{oo) on p and A. Color- 
coded values R{oo) from numerical simulations (a) and theoretical 
solutions (b) in the parameter plane (p, A), and the theoretical predic¬ 
tions are the solutions of Eqs. (|12^-jl4|( and l|16^-(|20^. Two vertical 
lines separate the plane into three regions, the former line is predicted 
by Eq. and the lat er line is predicted by numerically solving 
Eqs. l|21^, and j30[ ). In region I, only a vanishingly small frac¬ 
tion of individuals can be exposed to adopt the behavior (i.e., local 
behavior adoption). Region II shows a discontinuous phase transi¬ 
tion, while region III exhibits a continuous phase transition. The blue 
circles (blue dished solid curve) are the simulated (theoretical) crit¬ 
ical information transmission probability, respectively. We set other 
parameters as N = 10,000, p = 0.3 and Tb = 2. 


behavior above A^^, and further induce the bigots to adopt the 
behavior simultaneously above A^ [see black curves obtained 
from Eq. ([26]l]. In region IV (p > 0.5), half of the neighbors 
of bigots are activists. Once these activists adopt the behav¬ 
ior, the bigots will gradually adopt the behavior. Thus, R{oo) 
grows continuously and a finite fraction of individuals adopt 
the behavior above (see red curves). Our theoretical pre¬ 
dictions of A^ and R{oo) have a good agreement with the 
numerical predictions. 

It can be seen in Fig.|^that R(oo) increases discontinuously 
with A when Tf, = 2, thus R{oo) and phase transition on pa¬ 
rameter plane (p, A) for Tf, = 2 is explored in Fig.|^ And we 
find another crossover phenomenon in the phase transition; a 
change from being first-order to being second-order. Similar 
to Fig. I^is that the plane is divided into three regions: Re¬ 
gion I (p < 1 /9), the local behavior adoption region, in which 


only a vanishingly small fraction of individuals adopt the be¬ 
havior; region II (1/9 < p < 0.42) shows a first-order phase 
transition, where a finite fraction of individuals adopt the be¬ 
havior simultaneously above A^ (red dished lines); region III 
(p > 0.42) exhibits a second-order phase transition, in which 
R{oo) increases continuously versus A. The type of phase 
transition is verified by bifurcation analysis and studying vn. 


B. Heterogeneous networks 

We turn to elucidate the effects of network heterogeneity. 
To build the heterogeneous networks, the uncorrelated con¬ 
figuration model with power-law degree distributions P{k) ^ 
is adopted, where the mean degree {k) = 10 and the 
maximum degree k^ax ~ VN IWl . The network hetero¬ 
geneity increases with the decrease of 



















9 



0.0 0.2 0.4 0.6 0.8 1.0 

A 


FIG. 8 . (Color online) For Ti, = 4, the effects of network hetero¬ 
geneity on binary spreading threshold model, (a) R{oo) and (b) vr 
versus A on random regular networks, Erdos-Renyi (ER) networks, 
and scale-free networks with degree exponent 70 = 2.5 and 3.0. In 
figure (a), symbols represent the simulated results, and the lines are 
theoretical predictions from Eqs. l|12|>-jl4|l and (|16^-j20|>. Lines in 
(b) are the simulation results of vr. We set other parameters to be 
N = 10,000 andp = 0.3. 


In the case of Th = 4, we find that the global behavior 
adoption more likely to occur (i.e., lower in heteroge¬ 
neous networks, due to the existence of hubs in heterogeneous 
networks ll52l . as shown in Fig. Meanwhile, in strong het¬ 
erogeneous networks a large number of individuals with small 
degrees are difficult to adopt the behavior, so R{oo) is smaller 
at large A. For example, R{oo) at A = 0.7 is obversely smaller 
on scale-free networks with jd = 2.5 than that on RRNs. By 
bifurcation analysis of Eq. ( |2T| , it is discovered that the hybrid 
phase transition disappears for strong heterogeneous networks 
(e.g., 725 = 2.5 and 3.0 in Fig.|^. That is to say, network het¬ 
erogeneity leads to the crossover phenomenon: a change from 
being hybrid to being second-order. Moreover, the type of 
phase transition transition is further verified by observing V]i 
in Fig. ^b). Evidences in terms of the quantities R{oo), 
and A^^upport our edge-based compartmental theory. 

Eor the case of Tf, = 2 (see Eig.|^, we find the similar phe¬ 
nomenon of R{oo): R{oo) increases (decreases) with network 
heterogeneity for small (large) A. However, the system has the 
first-order phase transition, which denotes that network het¬ 
erogeneity does not alter the phase transition. Based on the 
bifurcation theory and the study of vr, the phase transition is 





FIG. 9. (Color online) For Tf, = 2, the effects of network hetero¬ 
geneity on binary spreading threshold model, (a) R(oo) and (b) vr 
versus A on random regular networks, Erdos-Renyi (ER) networks, 
and scale-free networks with degree exponent 70 = 2.5 and 3.0. In 
figure (a), symbols represent the simulated results, and the lines are 
the theoretical predictions from Eqs. and -l|20^. Lines 

in (b) are the simulation results of vr. And other parameters are set 
to be W = 10, 000 and p = 0.3. 


verified [see Eig.j^b)]. 


V. CONCLUSIONS 

Understanding social contagion dynamics in human popu¬ 
lations is extremely challenging. In practical behavior spread¬ 
ing, individuals usually display different criterions (wills) to 
adopt the behavior. That is to say, the heterogeneity of adop¬ 
tion thresholds do indeed exist, but its effects on social conta¬ 
gions have not been verified straightforward. To fill this gap, 
we proposed a non-Markovian behavior spreading model, in 
which individuals have distinct adoption thresholds, to ex¬ 
plore how heterogeneous adoption thresholds affect the final 
adoption size and phase transition. An edge-based compart¬ 
mental theory is developed to quantificationally describe this 
model, and this suggested theory is verified by a large num¬ 
ber of simulations. In the paper, we mainly focused on the 
so-called binary spreading threshold model, in which a frac¬ 
tion of individuals p have the adoption threshold = 1 and 
are acted as activists, and the remaining ones have a higher 
adoption threshold and are regarded as bigots. 
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We first studied the spreading dynamics on random regu¬ 
lar networks. And it is found that heterogeneous adoption 
thresholds markedly affect R{oo), and induce a hierarchical 
character in behavior adoption. In other words, activists first 
adopt the behavior and then stimulate bigots to adopt the be¬ 
havior. All the first-order, second-order and hybrid phase 
transitions are found to be in existence in the system. For 
the case of Tf, = 1, the traditional second-order phase tran¬ 
sition can be found. For Tf, > 3, the hybrid phase transi¬ 
tion mixed first and second order, i.e., R{oo) versus A first 
grows continuously and then follows by a discontinuous pat¬ 
tern, also occurs. More specifically, the continuous and dis¬ 
continuous growth of R{oo) are caused by the activists and 
bigots, respectively. There is a crossover phenomenon be¬ 
tween the hybrid and second-order phase transitions by vary¬ 
ing p. When Tfj = 2, the system only exhibits the first-order 
or second-order phase transitions. Interestingly, there is an¬ 
other crossover phenomenon; by varying p the phase transi¬ 
tion changes from being first-order to being second-order. 

Finally, we found that network heterogeneity markedly af¬ 
fects R{oo) and phase transition due to the existence of hub 
individuals. For the case of T{, > 3, strong network hetero¬ 
geneity causes the phase transition to change from being hy¬ 
brid to being second-order. For T}, = 2, network heterogene¬ 
ity does not alter the type of phase transition. 

The main contribution of our work lies in providing a quali¬ 


tative and quantitative view on the influence of heterogeneous 
adoption threshold. Meanwhile, our research results also en¬ 
rich the phase transition phenomenon. The developed theory 
can be generalized to behavior spreading with general adop¬ 
tion threshold distribution, and offer some new inspirations for 
other similar spreading dynamics, such as epidemic spread¬ 
ing and cascading. However, some fascinating and hopeful 
challenges still remain. For example, what will happen if the 
adoption thresholds are correlated with their degrees? How 
to extract more realistic behavior spreading mechanisms from 
real data? 
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